Mon. Not. R. Astron. Soc. OOP. [TlfT2l COOOO') Printed 11 July 2011 (MN I^TeX style file v2.2) 



The effects of thermohaline mixing on low-metallicity 
asymptotic giant branch stars 



o 
o 

(N 
o 

Q 

o 



(N 

o 

O 



X 



Richard J. Stancliffe* 

Centre for Stellar and Planetary Astrophysics, Monash University, VIC 3800, Australia 



Accepted 0000 December 00. Received 0000 December 00; in original form 0000 October 00 



ABSTRACT 

We examine the effects of tfiermolialine mixing on the composition of the envelopes of 
low-metallicity asymptotic giant branch (AGB) stars. We have evolved models of 1, 1.5 
and 2 from the pre-main sequence to the end of the thermally pulsing asymptotic 
giant branch with thermohaline mixing applied throughout the simulations. In agree- 
ment with other authors, we find that thermohaline mixing substantially reduces the 
abundance of '^He on the upper part of the red giant branch in our lowest mass model. 
However, the small amount of '^He that remains is enough to drive thermohaline mixing 
on the AGB. We find that thermohaline mixing is most efficient in the early thermal 
pulses and its efficiency drops from pulse to pulse. Nitrogen is not substantially af- 
fected by the process, but we do see substantial changes in -'^^C. The ^^C/^'^C ratio is 
substantially lowered during the early thermal pulses but the efficacy of the process is 
seen to diminish rapidly. As the process stops after a few pulses, the ^^G/^^C ratio is 
still able to reach values of 10'^ — 10^, which is inconsistent with the values measured 
in carbon-enhanced metal-poor stars. We also note a surprising increase in the ^Li 
abundance, with log]^Qe(^Li) reaching values of over 2.5 in the 1.5 Mq model. It is 
thus possible to get stars which are both C- and Li-rich at the same time. We compare 
our models to measurements of carbon and lithium in carbon-enhanced metal-poor 
stars which have not yet reached the giant branch. These models can simultaneously 
reproduced the observed C and Li abundances of carbon-enhanced metal-poor turn-off 
stars that are Li-rich, but the observed nitrogen abundances still cannot be matched. 
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1 INTRODUCTION 

It has long been known that models of asymptotic giant 
branch (AGB) stars that only include mixing in convec- 
tive regions are incomplete. These canonical models cannot 
account for observation s such as: the low ^^C/^^C ratios 
in low -mass AGB stars l|Abia fc Isern|[l997 : Lebzelter et al.l 
20081) . Li and C-ri ch stars in our Galaxy (|Abia fc Isernll 19971 
Uttenthaler et al] l2007h, isotopic r atios measured in pre- 
solax grains (e.g. NoUett et al.ll20()3l . and references therein). 
It has therefore been suggested that material might circulate 
below the base of the convective envelope and into regions 
where nuclear burning can happen. This process is often re- 
ferred to as 'cool bottom processing'. 

The carbon-enhanced metal-poor (CEMP) stars also 
show abundance trends that are difficult to explain in the 
context of AGB models including only convective mixing. In 
particular, those CEMP stars that are rich in s-process el- 
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ements (which we can, with some confidence, assume come 
from mass transfer in binary systems that once contained 
an AGB star) show ^'^C /^^C ratios no greater than 100 
(e.g. ISivarani et al.l '2006) , oven for stars at the Main Se- 
quence turn- off. AGB models pred ict ^'^C/^'^C ratios in ex- 
cess of 10" (IStancliffe et al.ll2009l ). In addition, the mod- 
els predict that low-mass AGB stars should be rich in car- 
bon, but nitrogen-poor. Yet there is a tendency for these 
stars to show nitrogen-enhancement alongside their carbon- 
enhancemen10. Furthermor e, there have been dete ctions of 
lithium in CEMP stars (e.g. [Thompson et aLlboOSl ). Canon- 
ical AGB models produce Li via the Cameron-Fowler mech- 
anism (Cameron & Fowler 1971), which involves the pro- 
duction of beryllium deep in the hydrogen burning shell via 
the reaction *He(^He,7)'^Be and the immediate transport of 



^ However. iMasseron et al] | |200g|) note that this trend only holds 
when considering CEMP stars as a whole. For each of the indi- 
vidual CEMP sub-classes, this trend does not hold. 
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this to cooler regions of the star where the ^hi that forms 
(once the beryllium has undergone electron capture) is sta- 
ble against proton captures. This takes place in the more 
massive stars which undergo hot bottom burning (HBB, 
where the base of the convective envelope lies in the top 
of the hydrogen-burning shell). Such stars would be rich in 
nitrogen, not carbon. These observations suggest that some- 
thing is missing from the AGB models and an extra mixing 
mechanism must be at work. 

Lithium is a particularly important element from 
the point of view o f mixing processes in CEMP stars. 
IStancliffe et al.l (|2007l ) pointed out that material accreted on 
to a low-mass companion does not just remain at the surface 
of the star and thermohaline mixing could efficiently mix 
this material deep into the stellar interior. However, detec- 
tion of lithium in the CE MP binary system CS 22964-161 led 
[Thompson et al] (|2008l ) to suggest that the mixing efficiency 
could not be so high. Li is a fragile element and is easily de- 
stroyed at temperatures in excess of about 2.5 x 10® K. Even 
a modest depth of mixing can lead to efficient Li-depletion 
l|Stancliffell200g| l and hence the measurement of Li in CEMP 
stars could be a good test of the efficiency of thermohaline 
mixing. It is therefore crucial that we understand the origin 
of this element. 

Despite the apparent need for extra mixing on both the 
giant branches, the physical nature of the mech anism (or 
mechanisms) has proved illusive. Recently, lEggleton et al.l 
showed that the lowering of the mean molecular 
weight by the reaction ^He(^He,2p)*He could lead to mix- 
ing in red giants via the thermohaline instability. This 
can potentially explain the change in abundances seen 
in gi a nts above the luminos ity bump jCharbonnel fc ZahnI 
l2007l : lEggleton et~aLl |2008| ). The effect of this mecha- 
nism has been investigate d beyond the first giant branch 
llCantiello fc Langeij lioog ) and in super-AGB stars l|Siessl 
I2009D . In this work, we wish to examine what the con- 
sequences of thermohaline mixing are for low-mass, low- 
metallicity AGB stars. 



2 THE STELLAR EVOLUTION CODE 

Calculations in this work have been carried out using the 
ST ARS stellar evolution code which was originally developed 
by lEggletonl (|l97ll 'l and has subsequ e ntly been updated by 
many authors (e.g. IPoIs et al] Il995l : IStancliffe fc Eldridgel 
l2009|). The v ersion used here includes the nucleosynthe - 



sis routines of IStancliffe et al.1 (|2005l ) and IStancUf^ ((2005I ). 
which follow the nucleosynthesis of 40 isotopes from D 
to ^■^S and important iron group elements. Th e code uses 



20041) ■ which em- 
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the opacity routines of lEldridge fc Tout! 
ploy interpolation in the OPAL tables 
1 19961 ) and which account for the variation in opacity as 
the C and O content of the material varies. An approxi- 
mation of the contribution to the molecular opacities is in- 
cluded via the method of Marigo ( 2002) and is described in 
IStanclifi'e fc Glebbeekl [20081^ 1. The modifications for foUow- 

^ Recently, low-temperature opacity tables for variabl e com- 
positions have becom e available l lLederer fc Aringej l2009t 
iMarieo fc Aringed I2OO9I) . These will eventually replace the less 
accurate approximations used here. 



ing the evolution through th e thermally-pulsing AGB (TP - 
AGB) phase are described in lStancliffe. Tout fc Poi3 \2004 ). 

Thermohaline mixing is included throughout 
all the evolutionary phases via t h e pre scription of 
iKippenhahn. Ruschenplatt fc ThomasI (Il980l ). with the 
mixing coefficient being multiplied by a fac t or of 100 as 
suggested by the work of lCharbonnel fc ZahnI (|2007t ). These 
authors find that with a factor of this magnitude they are 
able to reproduced the abundance t rends observed t o wards 
the tip of the red giant branch. IStancliffe et al.l (|2009l ) 
also showed that a coefficient of this magnitude could 
reproduced the observed mixing trends in both low-mass 
metal-poor stars and carbon-enhanced metal-poor stars on 
the upper part of the first giant branch. 

We evolve stars of 1, 1.5 and 2 M0 from the pre- main se- 
quence to the end of the thermally pulsing asy mptotic giant 
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branch (TP- AGB) using 999 mesh points. Reimerj (Il975l ) 
mass-loss prescription, with rj = 0.4, is used from the Main 
Sequence up to the TP- AGB; the lVassiliadis fc WoodI (|l993l ) 
mass-loss law is employed during the TP- AGB. A mixing 
length parameter of a = 2.0 is employed. The metallicity 
of each model is Z ^ 10^"* ([Fe/H]fh -2.3) and the ini- 
ti al abundances are assume d to be solar-scaled according 
to I Anders fc Grevessd (|l989l ). with the exception of ^Li, for 
which we adopt a value of Xtli = 1.05x 10~^ which is equiva- 
lent to the Spite plateau value. It would be more appropriate 
to adopt an a-enhanced composition for these models as the 
metal-poor stars are observed to have [a/Fe]« 0.4. However, 
opacity tables for an a-enhanced mixture that have also been 
computed with variable C and O abundances do not exist. 
We therefore choose to have our models self-consistent by 
keeping our initial compositions solar-scaled like our opac- 
ity tables. 



3 RESULTS 

The evolution of ^He, log^, e(''Li) and the ^^C/^^C ratio at 
the surface of each of the models, up to the beginning of 
the first thermal pulse in each model, is displayed in Fig. [1] 
The behaviour of each model is qualitatively similar. The 
''He abundance initially increases before dropping slightly 
and then levelling off. At the same time, the lithium abun- 
dance and the ^^C/^^C ratio both drop. These changes are 
all caused by the onset of first dredge-up - the deepening 
of the convective envelope as the star ascends the red giant 
branch. In each of the models, thermohaline mixing begins 
to affect the surface abundances at around 10** years before 
the first thermal pulse. The 1 Mq model is most strongly 
affected and its surface ^He abundance (by mass fraction) 
falls from around 9 x 10~* to about 1.5 x lO"**. At the same 
time, the remaining Li is efficiently destroyed by the mixing 
and the ^^C/^^C ratio drops from its post-dredge-up value 
of around 30 to around 5. 

The two higher mass models bring less ^He to the sur- 
face during first dredge-up. In addition, their surfaces are 
less affected by thermohaline mixing on the giant branch. 
There are only minor reductions in the surface abundances 



3 [A/B] = logio(AfA/AfB) - logio(AfA/AfB)Q, where A^, is the 
number abundance of species i. 
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Figure 1. Evolution of the surface abundances of various chem- 
ical species for the 1 Mq model (dotted line), 1.5 Mq model 
(dashed line) and 2M0 model (solid line). 



of ^He, ^Li and the ^'^C/^'^C ratio. The mean molecular 
weight inversion that formed during the giant branch is not 
completely removed by the time the star reaches the tip of 
the giant branch and consequently thermohaline mixing con- 
tinues to act throughout the core H e-burning phase (as was 
noted by ICantiello fc Langeijliooi ). The surface changes it 
causes are not dramatic and only the lithium shows a no- 
ticeable level of depletion. Eventually, as the convective en- 
velope deepens when the star ascends the early AGB, the 
surface abundances of ^He, ^Li and the ^^C/^^C ratio all fall 
slightly as the envelope once again dredges up material that 
has undergone nuclear burning. 

Each of the models is evolved to the onset of the super- 
wind phase at which point numerical problems occur and the 
runs were terminated. As it is clear that the effects we are 
interested in only occur at the beginning of the TP- AGB, 
we need not be concerned with this failure to complete the 
full thermally pulsing phase. In addition, because each of the 
models has entered the superwind, mass loss would quickly 
strip off the remaining envelope and we should not have 
missed more than a couple of pulse cycles. These missing 
pulses will not substantially affect the composition of the 
ejecta. Details of the models are given in Table [1] 

The occurrence of mixing by thermohaline convection 
driven by the burning of '^He is seen to occur in each of the 
models. Mixing is able to affect the abundances of ^He, ^Li 
and ^^C. Heavier elements are not affected as their proton- 
burning reactions are activated at much higher temperatures 
than the mixing mechanism is able to reach. We discuss in 
detail the case of the 1.5 M© model, the salient properties 
of which are displayed in Fig. (2] The other models behave 
in much the same way (see section 13.31 for a discussion of 
the differences between them). The mixing is at its most 
efficient in the earliest pulses. After the first thermal pulse 
with dredge-up, the ^^C/^'^C ratio falls dramatically in the 



interpulse. Just after the end of third dredge-up (TDUP) 
the ^^C/^^C ratio is around 600 while by the end of the 
interpulse it has fallen to just 120. At the same time the '^He 
abundances drops from 2.90 x 10"'' to 2.76 x 10~* and the ^Li 
abundances increases from zero to 1.99 x 10"^*^ (equivalent 
to logjg e(^Li) = l.SfJfl). On the next pulse, the effect is 
diminished with the ^^C/^'^C ratio being 540 at the end of 
TDUP and by the end of the interpulse it has dropped to 
350. The ^He abundances falls to 2.58 x lO""* and log^Q e(^Li) 
reaches 2.14. 

Why does the effect fall off so rapidly? The point at 
which the /x-minimum is located is critical. What deter- 
mines this is the competition between two things: the rate 
at which the mean molecular weight can be lowered via 
^He(^He,2p)'*He and the rate at which it is increased by 
other reactions producing "^He. In the absence of efficient 
CNO cycling in the H-shell the /i-minimum can lay deeper 
in the star as the pp-chains dominate. When dredge-up 
happens, the enhanced C-abundance increases the efficiency 
of the CNO cycle. This forces the /x-minimum outward to 
cooler temperatures and higher ^He abundances, with the 
magnitude of the /x-dip being reduced. The effect of this is 
two- fold. First, the reduction in the size of the /x-dip makes 
the mixing less efficient, as the /x-gradient (which the mix- 
ing rate is proportional to) is also reduced. As the location 
of the /i-dip occurs at progressively higher ^He abundances 
from pulse to pulse (the '^He abundance at the /x-minimum 
is around 6 x 10"^ during the second interpulse, 1.5 x 10~* 
during the third interpulse and 2 x 10"* during the fourth 
interpulse), the rate of ^He depletion falls off because mate- 
rial that is less depleted in ''He is now being mixed into the 
envelope. Secondly, the /i-minimum is pushed out to a re- 
gion of cooler temperature. In this region, the ^^C(p, 7)'^''N 
reaction is less efficient and the material that is mixed to the 
surface has a smaller fraction of ^"'C (formed from the decay 
of ^^N). This, coupled with the reduction in the mixing effi- 
ciency, means that the ^^C/^'^C ratio falls more slowly with 
each subsequent pulse. 



3.1 Lithium production 

Li production takes place in three stages and begins with 
the deepening of the convective envelope at TDUP. This ho- 
mogenises the convective envelope, resulting in a constant 
mean molecular weight above the soon to ignite H-burning 
shell. Without this flattening of the /x-profile, the burning of 
''He would not be able to produce a large enough change in 
the mean molecular weight to allow thermohaline mixing to 
take place. After TDUP has occurred, the convective enve- 
lope retreats, hydrogen burning reignites and the ^He-|-*He 
reaction leads to the production of ^Be (top panels of Fig.O. 
At the same time, the ''He(^He, 2p)'*He reaction leads to a 
reduction of the mean molecular weight. Once a /i- inversion 
develops, thermohaline mixing can take place and the beryl- 
lium is brought up to regions of cooler temperature where 
the ^Li that is formed from the beryllium is able to sur- 
vive for longer than it can closer to the burning shell (top 



logj^Q e(X) = logj^Q e(Afx/A'^H) + 12, where A^x is the number 
abundance of species X and Nn is the number abundance of hy- 
drogen. 
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Table 1. Details of the three models evolved. The columns are: TP - thermal pulse number; M« - total stellar mass in solar masses; 
Mh - hydrogen-exhausted core mass in solar masses; rip - interpulse period in 10^ years; logarithm of the maximum helium luminosity 
reached; AMh ^ growth of the H-exhausted core mass in solar masses; AMqup - mass of material dredged-up; A - dredge-up efficiency 
(which is defined as AMjjup/^Mh). The remaining columns give the abundances of carbon, lithium and the ^^C/^'^C ratio measured 
at the end of the interpulse phase following the thermal pulse. Each model run encountered numerical problems at the peak of a thermal 
pulse and hence no information for AMquPi a-i"! abundances could be given for the last pulse. 

right panel of Fig. [3)l . The profiles of ^Be and ^Li are due 
to the equilibrium between the rate at which these isotopes 
are transported up from the ''He burning regions (by the ac- 
tion of thermohaline mixing) and the rate at which they are 
destroyed. A pocket of ^Li is formed. As the /.t-inversion re- 
duces, the efficiency of mixing is reduced and ^Be is no longer 
mixed outwards as far. All this takes place very shortly after 
the cessation of TDUP. At no point is the efficiency of the 
mixing enough to get the Li into the convective envelope 
before it is destroyed. 

After the H-shell re-ignites the convective envelope 
moves back inward in mass (bottom left panel of Fig. [3]) 
and connects with the upper edge of the pocket of beryllium 
and lithium. This inward motion of the convective envelope 
takes place during the interpulse period and it is at this 
point that the Li enhancement is first seen. The envelope 



does not move in to as great a depth as it did during TDUP 
but it brings the base of the envelope close enough to the 
hydrogen-burning shell that Li can be transported up into 
the envelope before it is destroyed (bottom right panel of 
FigEJ. 



3.2 Carbon-13 

The production of ^''C is much more straightforward. Abun- 
dance profiles for the 1.5 model during its second inter- 
pulse phase are shown in Fig.U] with the isotopes ''He, '^^C, 
^■'C and ^''N being displayed. The newly-dredged up ^^C is 
converted into ^^C, which is then carried to the convective 
envelope by the thermohaline mixing. Once lifted from the 
H- burning shell, ^"'C undergoes no further reactions and so 
we see a steady increase in this isotope throughout the in- 
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Figure 3. Abundance profiles of the 1.5 Mq model after its second thermal pulse. In each panel, the abundances are: ''He (dot-short dash 
lino), '^Be (dot-long dash line) and ^Li (dashed line). The mean molecular weight, ^ is displayed by a solid line. The grey shading indicates 
convective regions. Top left: Just after the end of third dredge-up. Top right: Retreat of the convective envelope. Thermohaline mixing 
leads to the formation of a pocket of beryllium and lithium. Bottom left:. The hydrogen shell moves outward and the envelope comes 
back in again. Bottom right Just before the beginning of the next thermal pulse. 



terpulse (unlike ^Li, which still undergoes proton captures 
as it is being transported outward). Thus the ^^C/^^C ra- 
tio drops throughout the interpulse. There is a concomitant 
slight increase in the surface abundance as the temper- 
ature down to which mixing occurs is sufficiently high for 
incomplete CN-cycling. At the top of the H-burning shell, ni- 
trogen is enhanced to a few times its surface value after third 
dredge-up, whereas '^^C is over an order of magnitude more 
abundant at the top of the H-shell than it is at the surface. 
Hence the N-enhancement caused by thermohaline mixing 
is insignificant (and becomes even less significant with each 
pulse). 

During the third interpulse phase (see Fig. [5]), the 
12q^13q (^ggiijigg \yy rnuch less. The dredge-up of carbon 
forces the /x-minimum out to cooler temperatures. As a con- 
sequence of this, the CN cycle is less efficient and less ^"'C 
is produced: while the ^^C abundance in the envelope has 
increased from 3.1 x 10~* to 1.3 x 10~^, the ^^C abundance 
at the /i- minimum is roughly constant at around 3 x 10""'. 
In addition, the magnitude of the mixing coefficient is about 
a factor of 2 lower in the third interpulse than it was in the 
second, owing to the much smaller /x-dip. To compound all 



this, the interpulse is also 2 x 10* yrs shorter than the pre- 
ceding one. All this results in the ^^C/^^C ratio suffering 
less of a reduction during the interpulse. 

The reduction in the efficacy of mixing continues from 
pulse to pulse and very quickly thermohaline mixing is no 
longer able to affect the ^^C/^^C ratio. After three pulses 
there is barely a discernible change in it during the inter- 
pulse. The ^^C/^^C ratio keeps increasing with each episode 
of third dredge-up and values of over 10'^ are reached by the 
end of the TP- AGB. 



3.3 Mass variations 

The 2 M© model shows less of an effect on its surface abun- 
dances than the 1.5 M© model (see Fig.Ej). For the ^^C/^^C 
ratio, there is less time for circulation of material close to 
the H-burning shell because the interpulse periods are con- 
siderably shorter (by about a factor of 2) than in the 1.5 
model. In addition, the lower '^He abundance that this model 
starts the TP- AGB with also reduces the efficiency of mix- 
ing. The peak Li abundance that is reached is also lower 
than in the 1.5 M© model due to the reduced availability of 
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Figure 4. Abundance profiles of the 1.5 Mq model after its second tliermal pulse. In each panel, the abundances are: ''He (dot-short 
dash line), ^^C (dotted line), ^''C (short-dashed line) and (long-dashed line). The mean molecular weight, ^ is displayed by a solid 
line. The grey shading indicates convective regions. Top left: Just after the end of third dredge-up. Top right: Retreat of the convective 
envelope. Bottom left: The hydrogen shell moves outward and the envelope comes back in again. Bottom right Just before the 
beginning of the next thermal pulse. 



■'He and the reduced mixing efficiency. In addition, we also 
see a decrease in the Li abundance towards the end of the 
TP-AGB because the temperature at the base of the convec- 
tive envelope increases sufficiently to allow Li destruction. 
We also note that the Li abundance increases directly af- 
ter the first thermal pulse, even though there is no third 
dredge-up. This occurs because the convective envelope still 
moves inward after this thermal pulse but does not reach 
deep enough to penetrate the C-rich layers of the intershell. 
The envelope reaches to the H-burning shell and erases the 
existing mean molecular weight gradient. It is this flatten- 
ing of the mean molecular weight gradient by the convective 
envelope that allows the "'He burning to form a mean molec- 
ular weight inversion and thus allows thermohaline mixing 
to occur. 

The 1 Mq model also shows less of an effect on its sur- 
face abundances than the 1.5 M© model. Again, we attribute 
this to the lower ^He abundance (the model starts the TP- 
AGB with Xsjje = 1-39 X 10~*) which reduces the mixing 
efficiency. Unlike the 2 Mq model, the first efficient episode 
of thermohaline mixing happens after the first episode of 



third dredge-up. We thus get a substantial reduction of the 
i2Qyi3(2; j-atio during the next interpulse phase. 



4 UNCERTAINTIES 

There are many uncertainties associated with thermohaline 
mixing as a mechanism, most notable of which is the effi- 
ciency of the mixing. The prescriptions of lUlrichI (|l972l ) and 
iKippenhahn et all (|l980l ) give diffusion coefficients of the 
same functional form, but with a dimensionless constant (re- 
lated to the geometry of the mixing) that differs by orders of 
magnitude (see iCharbonncl fc Zahii 2007. for details). There 
are also differences in the way that one applies the mixing. 
lEggleton et al.l ||2008D use a mixing coefficient that depends 
on the difference between /i at a given point in the star 
and the mini mum value of n . This is a non-local prescrip- 
tion, whereas IStancliffe et~al] (I2007D use a local prescription 
where the mixing effi ciency depends on local differences in 
fj, (i.e. the u-gradient). lDenissenkov fc Pinsonneaultl l|2008bl ) 
also include overshooting in their prescription so that they 
actually mix to below the location of their /x-minimum. 
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Figure 5. Abundance profiles of the 1.5 Mq model after its second tliermal pulse. In each panel, the abundances are: ''He (dot-short 
dash line), ^^C (dotted line), ^''C (short-dashed line) and (long-dashed line). The mean molecular weight, ^ is displayed by a solid 
line. The grey shading indicates convective regions. Top left: Just after the end of third dredge-up. Top right: Retreat of the convective 
envelope. Bottom left: The hydrogen shell moves outward and the envelope comes back in again. Bottom right Just before the 
beginning of the next thermal pulse. 



In light of these uncertainties, we make the following 
tests on our 1.5 Mq model to determine the robustness of 
our results. We have run two model sets in which we change 
the diffusion coefficient of the existing prescription. In one, 
we increase the diffusion coefficient by a factor of fijfj (model 
D5) and in the other we decrease it by a factor of 10 (model 
DO.l). We also imple ment a prescription similar to that of 
lEggleton et~al] l|2008l) (hereafter EDL08), where the diffu- 
sion coefflcient is given by 

Ah = C(/i — ^imin) (1) 

for regions above the point of minimum fi. Here, Dth is the 
diffusion coefficient, C is a free parameter, fi is the mean 
molecular weight and /imin is the minimum mean molecular 
weight in the model. The free parameter C is set to 10^* as 
this gives approximately the same carbon depletion on the 
upper giant branch as is obtained in our standard model. Us- 
ing this prescription, the mixing is made global rather than 
local (model EDL). We note that the EDL08 prescription 

^ An increase of a factor of 10 proved to be unstable and the 
model could not be converged. 



gives a very different functional form for the mixing coeffi- 
cient to our standard implementation. The EDL08 prescrip- 
tion gives faster mixing further away from the /i-minimum 
because this is where n — Umin is greatest, whereas our stan- 
dard prescription gives most efficient mixing close to the 
point of /^-minimum because this is where the ^-gradient is 
at its steepest. 

Each of these models is run from the end of first dredge- 
up until the superwind phase of the TP- AGB. Table [2] shows 
how the ^Ife & ^Li surface abundances and the ^'^C/^^C ratio 
at various evolutionary stages change with the varying of the 
parameters. 

If we increase the diffusion coefficient by a factor of 5, 
we get substantially more processing of material on the first 
giant branch. The drop in the ^He during the mixing episode 
on the upper giant branch is about 6 times as great as for 
the standard diffusion coefficient. The lithium abundance 
drops from logj^g ^(^Li) = 0.41 to logj^g ^(^Li) = —0.3 and 
the ^^C/^^C ratio drops to just 10.1, compared to 21.2 in the 
standard case. Despite the depletion in the helium-3 reser- 
voir, the model is still able to undergo thermohaline mixing 
on the AGB. The total Li production is greater than in the 
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Figure 2. Properties of the 1.5 Mq model as a function of time 
since the first thermal pulse. Top panel: The evolution of the 
helium-burning luminosity. Middle panel: The evolution of the 
surface ^^C/^^C ratio (solid line) and the surface Li abundance 
(dashed line). Bottom panel: The evolution of the surface ^He 
abundance. 
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D5 


1.95 X 10"" 


-0.30 


10.1 


DO.l 


3.33 X 10-* 


0.35 


19.3 


EDL 


3.21 X 10"* 


0.32 


17.6 


First thermal pulse 


Standard 


2.76 X 10-4 


-1.40 


11.7 


D5 


1.45 X 10"* 


-1.82 


8.11 


DO.l 


3.03 X 10-* 


-0.48 


15.3 


EDL 


1.04 X 10-* 


-2.57 


6.02 


Final model 


Standard 


2.44 X IQ-* 


2.51 


2672 


D5 


1.22 X 10-* 


2.89 


2317 


DO.l 


2.92 X 10"* 


0.61 


10100 


EDL 


9.75 X 10-5 


0.24 


10420 



Table 2. Surface abundances by mass fraction of isotopes affected 
by thermohaline mixing and how they vary with parameters in 

the 1.5 Mq model, together with the ^•^C/"''^C ratio, at various 
evolutionary stages. The models are all started from the same 
point after the end of first dredge-up. 
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Figure 6. Properties of the 2 Mq model as a function of time 
since the first thermal pulse. Top panel: The evolution of the 
heliuin-buruiiig luminosity. Middle panel: The evolution of the 
surface ^^C/^^C ratio (solid line) and the surface Li abundance 
(dashed line). Bottom panel: The evolution of the surface ^He 
abundance. 



standard case by over a factor of 2, with the final model 
reaching a surface abundance of logj^Q e(^Li) = 2.89. The 
final ^^C/^^C ratio reached is 2317, which is only slightly 
lower than the standard case. This is because thermohaline 
mixing is only effective at reducing this ratio in the early 
pulses and this does not change with changes in the diffu- 
sion coefficient. 

Reducing the diffusion coefficient by an order of magni- 
tude does not substantially affect the amount of processing 
that takes place on the first giant branch. The surface abun- 
dances of •'He, ^Li and the *'^C/*''C ratio arc almost identical 
in the two models. The slower diffusion coefficient means 
there is less of a reduction in the ^He, ^Li and *^C/*^C 
ratio between the tip of the giant branch and the begin- 
ning of the TP-AGB. Along the TP-AGB, the reduction of 
the diffusion coefficient shows a unich more significant ef- 
fect. The reduction in the efficiency of mixing means that 
less '^Li is brought into the envelope and the model reaches 
only logj^Q e('^Li) = 0.61 by the onset of the superwind. The 
production of ^"'C is substantially reduced because of the 
weaker mixing and the model has a final ^^C/^'^C ratio that 
is a factor of about 4 higher than in the standard case. 

The EDL prescription gives similar changes in the abun- 
dances between the end of first dredge-up and the tip of the 
giant branch than those obtained in the standard case. This 
is not so surprising: we have chosen the free parameter of 
this prescription such that it gives roughly the same level of 
carbon depletion as the standard case, so we would expect 
the other species to behave in a similar way. However, the 
EDL prescription leads to significant depletion of ^He and 
'^Li at the beginning of the core helium burning phase. By 
the time the star reaches the AGB, the surface ^He abun- 
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dance has dropped to just 1.04 x 10~*, nearly a lactor of 3 
less than the standard case. This severely limits the action of 
thermohaline mixing along the AGB. With less ^He avail- 
able, the mixing efficiency is substantially reduced and so 
the drop in the ^He along the TP-AGB is much less than in 
the standard case. In addition, less ^^C can be transported 
up from the burning shell and the model reaches the super- 
wind phase with a very high ^^C/^^C ratio. ^Li production 
suffers two- fold. There is less ^He available from which to 
produce ^Be and the reduction in the transport efficiency 
means more beryllium and lithium is destroyed in the stel- 
lar interior before it can reach the safety of the convective 
envelope. The final surface lithium abundance is over 2 dex 
lower than in the standard case, with logj^Q e(^Li) being just 
0.28. 

There are also uncertainties that we are unable to as- 
sess at present. Foremost among these is the suggestion that 
thermo haline mixing can be inhibited by th e presence of ro- 
tation ijPenissenkov fc Pinsonneaultll2008al ) . Until hydrody- 
namic simulations of the interactions of these two mecha- 
nisms are made, this point cannot be addressed. It is also 
probable that other effects (e.g. the presence of magnetic 
fields) will also infiuence the degree of mixing. 



5 DISCUSSION 

The action of thermohaline mixing during the TP-AGB does 
reduce the ^'^C/^'^C ratio, so that the ejecta will have a 
ratio of a few thousand rather than the 10^ predicted by 
canonical models (e.g. Stancliffc & Globbook 2008). While 
this is an improvement it still leaves the models predicting 
higher ^^C/^^C ratios than are observed in CEMP turn-off 
stars. While more evolved CEMP stars will undergo mix- 
ing of CN-cycled material on the first giant branch (either 
at first dredge-up or through extra mixing above the lumi- 
nosity bump) which will lower their surface ^^C/^^C ratios, 
the turn-off objects are not evolved enough for this to have 
happened. Another mechanism for reducing this ratio must 
clearly be sought. It is also unable to affect the abun- 
dance and so cannot account for the correlation of C- and 
N-enhancement observed in the CEMP starfl However, if 
mixing occurs to regions belo w the point of minimum due 
to som e sort of overshooting (|Denissenkov fc Pinsonneaultl 
l2008bl ). it is possible that substantial enhancements of ^'^C 
and could result. 

The surprising outcome of these simulations is the high 
Li abundances that can be produced. It is usually supposed 
that only the higher mass AGB stars which undergo hot bot- 
tom burning are able to produce Li v ia the Cameron- Fowler 
mechanism (iCameron fc Fowleij|l97ll ) . This work shows that 
it is possible that low-mass AGB stars could be producers 
of lithium-7. We also note that the action of thermohaline 
mixing on the AGB is subtly different from its action of the 
RGB. On the RGB, it leads to a depletion of Li with the 
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[Fe/H] 


[C/Fe] 


logio "^CLi) 


log g 


Refs. 


CS 22964-161A 


-2.41 


1.35 


2.09 


3.7 


1 


CS 22964-161B 


-2.39 


1.15 


2.09 


4.1 


1 


HE 0024-2523 


-2.7 


2.6 


1.5 


4.3 


2 


CS 31080-095 


-2.85 


2.69 


1.73 


4.5 


3 


CS 31062-012 


-2.53 


2.14 


2.3 


4.3 


4 



Table 3. Properties of CEMP turn-off stars with measured Li- 
abu ndances, as extrac t ed fro m the SAGA database . Refer ences: 
1 - iThompson et al.l | |2008|). 2 - iLucatello et al.l (l2003f) . 3 - 



ISivarani et al.l ||2006|) . 4 - lAoki et al.l ||2008|') " 



star leaving the RGB with virtually no lithium left. How- 
ever, on the AGB thermohaline mixing can substantially in- 
crease the Li abundance up to values above the Spite plateau 
value. It is therefore possible that low-mass AGB stars have 
contributed to the Galaxy's Li budget. The effect of a pop- 
ulation of low-mass, low-metallicity lithium producers on 
Galactic chemical evolution models should be investigated. 

These models do improve the agreement of the AGB 
models with observations of Li in carbon-enhanced metal- 
poor stars. We have extracted from the Stellar Abundances 
for G alactic Archaeology (SAGA) database (|Suda et al.l 
l2008f) those CEMP stars that are both C-rich and have 
measured Li-abundances, and also that are still close to the 
main sequence turn-off (because first dredge-up will signif- 
icantly reduce the surface Li abundance). We select only 
those stars in the metallicity range —3 <[Fe/H]< —2 as 
the models presented herein may be expected to apply only 
over a limited range in metallicity. In particular, at very 
low metallicities, AGB stars can undergo ad ditional mixing 
events not present in our models (see e.g. Fuiimoto et al.l 
Il99d : ICampbeU fc Lattanzid l2008l : iLau et aLlBo09l . among 
many others for a discussion of these mixing events). SAGA 
lists 5 turn-off objects with measured lithium abundances, 
the properties of which are displayed in Table |31 

One caveat should be added to the following discussion. 
The scenario of mass transfer from an AGB primary star on 
to a lower mass secondary in a binary system is expected to 
apply to those stars belonging to the CEMP-s subclass, i.e. 
those stars which have [Ba/Fe]> 1 and [Ba/Eu]>0.5 accord - 
ing to the definitions given by iBeers fc Christliebl (|2005h . 
The origin of the r/s subclass of CEMP stars, which have 
<[Ba/Eu]< 0.5, is currently unknown. It is possible that 
their s-process enrichment has come from a binary mass 
transfer event, in which case the models presented herein 
would apply, b ut the enrichm ent may have another source 
entirely (see e.g. lLugaro et alj 2009. for a possible alternative 
formation scenario). This should be borne in mind through- 
out the following discussion. 



^ iMasseron et al.l | |2009| ) suggest that the apparent C and N cor- 
relation does not hold when considering the individual CEMP 
subclasses. Therefore, the failure of the model to produce a C 
and N correlation may not be so serious. 



5.1 CS 22964-161A & B 

iThompson et al] l|2008l ) reported abundances for this in- 
triguing system. It is a double CEMP binary, with both 
components being turn-off objects. It is presumed that this 
system is actually a hierarchical triple system, in which an 
AGB star polluted the tight inner binary which we now ob- 
serve today. 

As both objects in the system have [C/Fe] = 1.21 we are 
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forced to assume that some mixing of accreted material has 
taken place in order to get the surface [C/Fe] low enougfQ, as 
the three A GB mo dels presented here all give [C/Fe]> 2 in 
their electa. [Stancliffe ( 2009.) looked at models in which ther- 
mohaline mixing and gravitational settling are allowed to 
occur. In additio n, an extra turbulent m ixing process might 
also be at work. iRichard et al] (|2005l ) suggested that the 
Spite plateau value for Li could be reconciled with the Big 
Bang Nucleosynthesis predictions if some unknown turbu- 
lent mixing process near the stellar surface could carry Li 
down to temperatures where some of it could be destroyed. 
The physical cause of such a process remains elusive and the 
prescription employed is an ad hoc one. 

Models of the se condary involving these three processes 
fsee lStancliffil2009l . for details) were run, accreting 0.001- 
0.1 Mq of material from companions of 1, 1.5 and 2Mq, 
in each case producing a star of 0.8 Mq. The evolution of 
these models through the [C/Fe]-logjg e(Li) plane are shown 
in the top panel of Fig. [T] We find that we can reproduce 
the observed properties of OS 22964-161 if around 0.002 Mq 
of material was accreted from the 1.5 Mq companion, with 
star A accreting slightly more mass than star B. In contrast, 
models including thermohaline mixing alone deplete far too 
much of their lithium to ma tch the ob s ervati ons (lower panel 
of Fig. 0, as was noted bv lStancliff j l|2009l ). 



5.2 HE 0024-2523 

Detailed abun dance analyses for th is ob ject have been 
carrie d out by ICarretta et al.l (|2002l ) and iLucatello et al.l 
l|2003l '). These authors report [Fe/H] for this object as be- 
ing -2.62 and -2.72 respectively. Both studies agree that the 
surface gravity of the object is logg = 4.3. The bulk of the 
abundances for this object come from the latter work and 
it is from here that we quote the abundances [C/Fe] =2. 6 
and login efLi) = 1-5. The barium abundan ce according to 
ILucatell o et al.l l|2003l ) is [Ba/Fe]=1.46, while lCarretta et all 
(2002) give [Ba/Fe] = 1.76. The object is therefore s-process 
rich. Onl y an upper lim i t for europium, [Eu/Fe]< 1.16 is 
given by ILucatello et al.l l|2003l ) and so it cannot be ruled 
out that HE 0024-2523 belongs to the r/s CEMP sub-class. 
However, these authors do report that the star is extremely 
lead-rich, with [Pb/Fe] = 3.3. 

HE 0024-2523 could conceivably have accreted material 
from a companion in a similar mass range to the models pre- 
sented here. The [C/Fe] and logj^g e(Li) values of the models 
are in the region of the values observed in this object. We 
would have to assume that no mixing of the accreted mate- 
rial took place, as any mixing would reduce the C- and Li 
abundances too much. 



5.3 CS 31080-095 

ISivarani et al.l (|2006t ) have performed a detailed abundance 
analysis of this object. They find [Fe/H] = -2.85, [C/Fe] = 



^ Alternatively, the AGB models could be over-predicting the 
amount of C that is dredged-up. A simple back of the envelope 
calculation suggests that to reach [C/Fe] = -1-1 about 5 X 10"'' Mq 
of intershell material would have to be dredged-up, compared with 
around 4 X 10~^ Mq in the 1 Mq model. 
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Figure 7. The evolution of logjQ e(Li) with [C/Fo] when accret- 
ing material from a 1.5 Mq companion. The cases displayed are 
for when 0.001 Mq (solid line), 0.01 Mq (dotted line) and 0.1 Mq 
(dashed line) is accreted. In each case, the secondary is left with 
a total mass of 0.8 Mq. Bold lines indicate where logg passes 
from 4.5 to 3.5 as the object evolves off the main sequence. The 
errorbars denote the locations of specific observed systems. Top 
panel: The secondary is modelled including thermohaline mix- 
ing, gravitational settling and an extra turbulent process. Bot- 
tom panel: The secondary is modelled with thermohaline mixing 
alone. 



2.69 and logj^Q e(Li) = 1.73. A moderate barium enhance- 
ment of [Ba/Fe] = 0.77 is reported, but this is the heaviest 
element for which they give an abundance. 

CS 31080-095 could conceivably come from a similar 
mass range to the models presented here. The C and Li 
abundances of this object lie between those of our 1 and 
1.5 Mq models, so we would expect that accreting material 
from a companion in this mass range would give appropriate 
abundances, provided the accreted material did not mix into 
the secondary. 

5.4 CS 31062-012 



CS 31062- 12 ha s been stud i ed by iNorris et al] lll997h 
lAoki et all (l200ll ) lAoki et al.l l|2002al ). lAoki et al.l l|2002bD 
and lAoki et al.l (120081 ). These authors give [Fe/H] in 

the range - 2.53 to -2 .74, with [C/Fe] in the range 

2.1 to 2.15. lAoki et all (120081 ) give the lithium abun- 
da nce as logi n e(Li) = 2.3, but an earlier study 
by ICharbonnel fc Primad (|2005l ) gives a lower value of 
logjQ e(Li) = 1.973. The object is barium-rich, with [Ba/Fe] 
around 2, and somewhat enriched in europium w ith [E u/Fe] 
= 1.4 bei n g re ported by both iNorris et al.l ^9M) 

and 

lAoki et all (|2002b). This would place CS 31062-012 in the 
CEMP-s subclass. 

We find we are able to model this system using ejecta 
from a 1.5 Mq companion. If around 0.002 Mq of material is 
accreted on to the secondary, and both thermohaline mix- 



Thermohaline mixing in low-Z AGB stars 11 



Object 


[Fe/H] 


[C/Fe] 


logio "^CLi) 


log g 


Refs. 
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2.72 
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2.1 
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4.2 


2 


SDSS 2047+00 


-2.05 


2.3 


< 2.3 


4.5 


2 



of the accreted material is a viable scenario. For CS 29497- 
030, accretion from a 1 M© companion is consistent with the 
observed abundances, provided no mixing of the accreted 
material takes place. However, it is also possible the com- 
panion was more massive and that mixing of the accreted 
material happened. 



Table 4. Properties of CEMP turn-off stars with upper limits on 
the Li ab undance, as ext r acted fro m the SAGA database. Refer- 
ences: 1 - ISivarani et al.l l l2004h . 2 - lAoki et all ||2008|) . 

ing and gravitational settling are included (but not the ad- 
ditional turbulent process), then the properties of the sys- 
tem can be reproduced. Dilution of the small quantity of 
material accreted allows the carbon abundance to drop to 
the requisite level, while the gravitational settling creates 
a mean molecular weight barrier that prevents the mixing 
going too deep, allowing a substantial quantity of lithium to 
survive. Inclus i on of the extra turbulent mixing process of 
iRichard et aD (|2005l ) leads to a degree of carbon depletion 
that is too great. 

It is rather unsatisfactory that for some objects we seem 
to require some degree of mixing of the accreted material 
and for others we require that no mixing of thi s material 
taJces place. Similar conclusions were reached bv IStancliffd 
when examining the trends of various light elements 
in CEMP stars. The physical reason why some stars mix 
their accreted material and others do not remains unknown. 
It is likely that the interaction of these mixing processes 
with other physical processes not included here (e.g. rota- 
tion, magnetic fields) is responsible. 

5.5 Systems with upper limits on lithium 

In addition to the above systems which have measurements 
of the lithium abundances, SAGA lists five systems for which 
an upper limit for lithium has been determined. These sys- 
tems are: CS 29497-030, CS 29526-110, SDSS 0924+40, 
SDSS 1707+58 and SDSS 2047+00. The salient details of 
these models are given in Table [J] The high upper limits 
derived for the latter four systems do not severely restrict 
which models can be applied. For example, a model involving 
accretion of material from a 1 M0 AGB star, with no mixing 
of the accreted material into the secondary would give sur- 
face abundances broadly consistent with the observations of 
all but SDSS 0924+40. However, accretion of ejecta from a 
more massive AGB companion would also fit the data pro- 
viding some mixing of the accreted material took place in or- 
der to lower the surface carbon abundance of the secondary. 
Without better constrains on the lithium abundance, it is 
not possible to determine which of these scenarios actually 
took place. 

The systems CS 29497-030 and SDSS 0924+40 are more 
restrictive. The high [C/Fe] of SDSS 0924+40 rules out ac- 
cretion from a 1 M0 companion because this model does 
not produce enough carbon. We can also rule out the case 
of accretion of material from a 1.5 M© model. While this 
model produces roughly the correct amount of [C/Fe], its Li 
abundance is about 0.5 dex above the observed upper limit. 
Accretion from a more massive companion with some mixing 



6 CONCLUSIONS 

We have investigated the effect that thermohaline mixing 
has on the abundances of low-mass, low-metallicity AGB 
stars. We find that enough "^He remains after the first giant 
branch that thermohaline mixing can still take place on the 
AGB. However, the effect is only felt during the very first 
thermal pulses with dredge-up. Thermohaline mixing can 
lead to substantial production of ^Li - even up to values 
above the Spite plateau value. Thus it is possible to recon- 
cile C- and Li-rich metal-poor stars with having come from a 
binary mass transfer scenario. We demonstrate that the Li- 
enrichments of five known Li-rich, turn-off CEMP stars can 
be explained using the abundances of our models. Thermo- 
haline mixing does not reduce the ^^C/^'^C ratio sufficiently 
to reconcile the AGB models with observations of CEMP 
stars at the Main Sequence turn-off, nor does it raise the 
surface nitrogen abundance sufficiently. The possibility that 
low-mass, low-metallicity stars could be producers of lithium 
is intriguing and their role in Galactic chemical evolution 
should be investigated. 

Many uncertainties about the nature of thermohaline 
mixing still remain and we have shown that the Li abun- 
dance is particularly sensitive to the nature of the mixing. 
There is also the question of how other physical processes 
(e.g. rotation) interact with thermohaline mixing and the 
consequences of this remain to be explored. 
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